Asymptotic estimation theory 
for a finite dimensional pure state model 



Masahito Hayashi 
Department of Mathematics, Kyoto University 

Kyoto 606-8502, Japan 
e-mail address: masahito@kusm.kyoto-u.ac.jp 

Abstract 

The optimization of measurement for n samples of pure sates are studied. The error of the 
optimal measurement for n samples is asymptotically compared with the one of the maximum 
likelihood estimators from n data given by the optimal measurement for one sample. 



1 Introduction 

Recently, there has been a rise in the necessity for studies about statistical estimation for the 
unknown state, related to the corresponding advance in measuring technologies in quantum optics. 
An investigation including both quantum theory and mathematical statistics is necessary for an 
essential understanding of quantum theory because it has statistical aspects [Q, ||. Therefore, 
it is indeed important to optimize the measuring process with respect to the estimation of the 
unknown state. Such research is known as quantum estimation, and was initiated by Helstrom in 
the late 1960s, originating in the optimization of the detecting process in optical communications 
[0 . In classical statistical estimation, one searches the most suitable estimator for which probability 
measure describes the objective probabilistic phenomenon. In quantum estimation, one searches 
the most suitable measurement for which density operator describes the objective quantum state. 

Contained among important results are three estimation problems. The first is of the complex 
amplitude of coherent light in thermal noise and the second is of the expectation parameters of 
quantum Gaussian state. The former was studied by Yuen and Lax j^] and the latter by Holevo 
||. These studies discovered that heterodyning is the most suitable for the estimation of the 
complex amplitude of coherent light in thermal noise. The third is a formulation of the covariant 
measurement with respect to an action of a group. It was studied by Holevo || Q. In the 
formulation, he established a quantum analogue of Hunt-Stein theorem. 

Quantum estimation, was first used in the evaluation of the estimation error of a single sample 
of the unknown state as it had advanced in connection with the optimization of the measuring 
process in optical communications. Thus early studies were lacking in asymptotic aspects, i.e. 
there were few researches with respect to reducing the estimation error by quantum correlations 
between samples. 

Recently, studies about the estimation of the unknown state are attracting many physicists 
[§i @; Hi ■ Some of them were drawn by the variation of the measuring precision with respect to 
the number of samples of the unknown state ||, . 

Nagaoka [p"I| studied, for the first time, asymptotic aspects of quantum estimation. He paid 
particular attention to the quantum correlations between samples of the unknown state, and studied 
the relation between the asymptotic estimation and the local detection of a one-parameter family 
of quantum states. 

In the early 1990s, Fujiwara and Nagaoka [l^, [u] studied the estimation problem for a multi- 
parameter family consisting of pure states. They pioneered studies into the estimation problem 
of the complex amplitude of noiseless coherent light. The research found that heterodyning is the 
most suitable for the estimation of the complex amplitude of noiseless coherent light as for the one 
of coherent light in thermal noise. In 1996, Matsumoto jl5) established a more general formulation 
of the estimation for a multi-parameter family consisting of pure states. Moreover in 1991, Nagaoka 
J16| treated the estimation problem for 2-parameter families of mixed states in spin 1/2 system, 



and in 1997 the author jl7j [l8| treated it for 3-parameter families of mixed states in spin 1/2 
system. However, there are no asymptotic aspects in these works about multi-parameter families. 
There is more necessity of this type of investigation into one- and multi-parameter families. 

Can quantum estimation reduce the estimation error by using the quantum correlations between 
samples, under the preparation of sufficient samples of the unknown state? To answer this question, 
in this paper, we treat a family, consisting of all of pure states on a Hilbert space TL [] under the 
preparation of n samples of the unknown state, with the estimation problem In §|^, we use, as a 
tool, the composite system consisting of n samples as a single system. The quantum i.i.d. condition 
is introduced as the quantum counterpart of the independent and identical distributions condition 
(|3|). In §^|, we review Holevo's result concerning covariant measurements which will be used in the 
following sections. In §^|, we apply Holevo's result to the optimization of measurements on the 
composite system, which results in obtaining the most suitable measurement (Theorem ||). We 
asymptotically calculate the estimation error by the optimal measurement in the sense of both the 
mean square error (MSE) and large deviation, (see (PI) (|lC|) (|ll|) (|l3|).) The first term of the right- 
hand side of (|l^) is consistent with the value conjectured from the results by Fujiwara, Nagaoka 
E3 and Matsumoto Jl5| ). However, the optimal measurement may be too difficult for modern 
technology to realize when using more than one samples. 

In §[|, we use this estimation problem under the following guidelines. The samples are divided 
into pairs consisting of a maximum of m samples. By measuring each pair with the optimal mea- 
surement in section ^, we create some data. The estimated valued is given by manipulating these 
data. The restricted condition is called m-semiclassical (see (|l4|)). We compare an m-semiclassical 
measurement with the optimal measurement of section |l|with respect to the estimation error under 
the preparation of a sufficient amount of samples. When we use the maximum likelihood estimator 
to manipulate the data, the MSE of both asymptotically coincide in the first order (see (|l0|)(^9|)). 
However, when the radius of allowable errors is finite, the error of large deviation type in the latter 
is smaller than that in the former type (see (pd])(pO|)). 

Can we asymptotically realize a small estimation error as the optimal measurements in section \\ 
has? It is, physically, sufficient to construct the optimal measurement for one sample. In section |g, 
we show how to construct it (see (p5|)). And in this case, we can calculate the maximum likelihood 
estimator from data by using computer. 

Most of the proofs of this paper are given in Appendices. In view of multiparameter families of 
mixed states in spin 1/2 system, Hayashi has discussed the same problem by using Cramer- Rao 
type bound. 



2 Pure state n-i.i.d. model 

In this section, we use the mathematical formulation of the estimation for pure states. Let k be 
the dimension of the Hilbert space TL, and V(TL) be the set of pure states on TL. 

In quantum physics, the most general description of a quantum measurement is given by the 
mathematical concept of a positive operator valued measure (POVM) D, I on the system of state 
space. Generally, if f2 is measurable space, a measurement M satisfies the following: 

M (B) = M(B)*,M(B) > 0, M(0) = 0, M(O) = Id on TL, for any B C fl. 
M(UiBi) = M ( B i)> for B i n B i = ^ i): i B i\ is countable subsets of fi. 

i 

In this paper, M. (J7, TL) denotes the set of POVMs on TL whose measurable set is fi. A measurement 
M £ M.(Q, TL) is called simple if M{B) is a projection for any Borel B C f2. A measurement M is 
random if it is described as a convex combination of simple measurements. A random measurement 
M = ^ i aiMi {Mi is simple and > 0.) can be realized when every measurement Mi is done 
with the probability a^. 

1 Where H denotes a finite-dimensional Hilbert space which corresponds to the physical system of interest. 



2 



In this paper, we consider measurements whose measurable set is V(H) since it is known that 
the unknown state is included in V(H). 

Next, we define two distances charactering the homogeneous space V(H). 

Definition 1 the Fubini-Study distance df s (which is the geodesic distance of the Fubini-Study 
metric) is defined as: 

cosd fs (p,p) = ^tvpp , < d fs (p,p) < -. (1) 
the Bures distance d\, is defined in the usual way: 



d b (p,p) ■= tr pp. (2) 

It is introduced by Bures \2f\J in a mathematical context. 



Let W(p,p) be a measure of deviation of the measured value p from the actual value /?, then we 
have the following equivalent conditions: 

o W(p,p) = W(gpg*,gpg*) for g e SV(k),p,p e V(H). 

o There exists a function h on [0, 1] such that W(p, p) = ho df s (p, p). 

It is natural to assume that a deviation measure W(p,p) is monotone increasing with respect to 
the Fubini-Study distance df s . 

If 7i\, . . . 7 H n are n Hilbert spaces which correspond to the physical systems, then their com- 
posite system is represented by the tensor Hilbert space: 

U {n) ;=Hi® ■■■®H n = ® Hi. 

i=l 

Thus, a state on the composite system is denoted by a density operator p on H^ . In particular if 
n element systems {Tti} of the composite system TiS n ' are independent of each other, there exists 
a density pi on Hi such that 

, , n 
p (n> = pi <8> • • • ® p n = <8> Pi. 

i=l 

The condition: 

Hi = ■ ■ ■ = H„ = H, pi = ■ ■ ■ = p n = p (3) 

corresponds to the independent and identically distributed condition (i.i.d. condition) in the clas- 
sical case. In this paper, we treat with this estimation problem under this condition (j^) called 
the quantum i.i.d. condition. This condition means that identical n samples are independently 
prepared. The model {p^ n ' — p® ■ ■ ■ ® p\p G V(H)} is called rt-i.i.d. model. As p is a pure state, 

n 

H {n) and p( n ) are simplified as follows. Letting p = \</>){<f>\ € V(H), we have 



Because all of the vectors (j)^ is included in n-times symmetric tensor space, for any mea- 
surement M € Ai(fl, H^) on the n-times tensor space 7v- n ', the measurement M{du>) := 
P^An) M ( du) P^An) € A4(Q, Hi ) on the n-times symmetric tensor space Hs satisfies that: 

txM(dw)p^ = trM(duj)p {n) for any peH, 

(n) 

where H s denotes the n-times symmetric tensor space on H. Therefore, all of possible mea- 
surements can be regarded as elements of M(V(H), Hi ). The mean error of the measurement 
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II G M(V(H),H S " ) with respect to a deviation measure W(p,p), provided that the actual state 
is p, is equal to 

Pf'W(n):= / ^(p,p)tr(n(rfp)pW). 

In minimax approach the maximum possible error with respect to a deviation measure W(p, p) 

X>w,(i»)m) : = max Df'Wfll) 

is minimized. 

3 Quantum Hunt- Stein theorem 

In this section, the quantum Hunt-Stein theorem, established by Holevo 0), is summarized. Let 
G be a compact transitive Lie group of all transformations on a compact parametric set 0, and 
{V g } a continuous unitary irreducible representation of G in a finite-dimensional Hilbert space 
W := C fe ', and [i a ct- finite invariant measure on group G such that p(G) = 1. In this section, we 
consider the following condition for a measurement. 

Definition 2 A measurement II G A4(Q,7i') is covariant with respect to {V g } if 

v;n(B)v^ = n(vO 

/or any jeG and any Borel Bc9, where 

B g := {g9\9 G £?}. 

.M(0, V) denotes the set of covariant measurements with respect to {V g }. 
Covariant measurements are characterized by the following theorem. 

Theorem 1 The map V s from the set S(7i') of densities on TV to Ai(Q, V) is surjective for any 
9 € G, where V (P) is defined as follows: 

V d (P)(B) := k' f VgPV g *p(dg) for B G B(Q) 
J{ g eeB} 

for any P £ S{W). 

In this section, we treat with the following condition for a family of states. 

Definition 3 The family is called covariant under the representation {V g } of group G acting on 
&,if 

S g e = VgS e V; i V 3 eG,V(Je9. 

Assuming that the object is prepared in one of the states {Sg\6 € 0} but the actual value of 9 is 
unknown, then the difficulty is estimating this value as close as possible to a measurement on the 
object. We shall solve this problem by means of the quantum statistical decision theory. 

Let W(9, 9) be a measure of deviation of the measured value 9 from the actual value 9. It is 
natural to assume that W(9,9) is invariant: 

W{9, 9) = W(g9, g§) for Vg G G, V0, V(? G 0. (4) 

The mean error of the measurement II G Ai(Q,H') with respect to a deviation measure W(9, 9), 
provided that the actual state is Sg, is equal to 

£>f' 5 (IT) := f W(6,6)tr(IL(d§)Sg). 
Je 
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Following the classical statistical decision theory, we can form two functionals of T>$ giving a total 
measure of precision of the measurement II. 

In Bayes' approach we take the mean of T>^ with respect to a given prior distribution ir( d9). 
The measurement minimizing the resulting functional: 

2?f< s (n) := / vj' s {xi)n{de) 

Je 

is called Bayesian. This quantity represents the mean error in the situation where 9 is a random 
parameter with known distribution n(d9). In particular, as 0,G are compact and "nothing is 
known" about 6, it is natural to take for ir( d6) the "uniform" distribution, i.e. normalized invariant 
measure v{ d9) defined as follows: 

u(B) := n{{ 9 6 E B}). 

It is independent of the choice of 9 S G. 

In minimax approach the maximum possible error with respect to a deviation measure W(9, 9) 

V W > S (I1) :=max2?^ S (n) 
y J see e K 1 

is minimized. The minimizing measurement is called minimax. 

Because G is compact, we shall show that in the covariant case the minima of Bayes and 
minimax criteria coincide and are achieved on a covariant measurement. We obtain the following 
quantum Hunt-Stein theorem [|[ It is easy to prove the theorem. 

Theorem 2 For a covariant measurement II G M(0, V), we obtain the following equations: 
For ne M(Q,W), denote 

U g (B) := V g n(B g )V g * for B € B(Q). 
Introducing the "averaged" measurement 

fl(B) := f VW*). 
Jg 

we have 

2^ s (n)= [ v^ s (n g - 1 Md 9 ) = v^ s (n). 

Jg 

Thus, 

V w ' s {n) > £>r< s (n) = V^ S (U). 

In this case, minimax approach and Bayes' approach with respect to v(d9) are equivalent. There- 
fore we minimize the following: 

Vj' s oV e (P) = k' f W(0 >9 e)tTS e V g PV g *n(dg) = tvW(e)P, 
Jg 

where 

W{9) := k' f W{9,g9)V;S e V g ^{dg) 
Jg 

= k' [ W(6,6)S § v{d§). 
Je 

Thus, it is sufficient to consider the following minimization: 

min tvW(e)P= min txW(6)P. 
PeS(H) Pev(H') 
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4 Optimal measurement in pure state n-i.i.d. model 

In this section we apply the theory of ^ to the problem §0. 
We let as follows: 

9 := V(H), H' := Tif, G := SU(fc), S p := 

Then, the invariant measure v on V(H) is equivalent to the measure defined by the volume bundle 
induced by the Fubini-Study metric. We let the action {V g } of G = SU(fc) to Hg be the tensor 
representation of the natural representation. In this case, we have k' = 

Theorem 3 If a deviation measure W(p, p) is monotone increasing with respect to the Fubini- 
Study distance df s , then we get 

min tr W(p)P = tr W{p)p [n) . 

For a proof see Appendix [A]. Thus, V p (p^) is the optimal measurement with respect to a devi- 
ation measure W(p,p). The optimal measurement is independent of the choice of p and W since 
V po (pQ ) = V p (p( n ^). This optimal measurement is denoted by H n and is described as follows: 



n„(dp) 



n + k-1 
k-l 



p^u(dp). 



Under the following chart (||), the optimal measurements are denoted as: 



n n (de) = 



n + k — 1 
k-l 



(t>{o) {n) ) U(ei n) v(do) 



(•5) 



for 9 £ {9 £ R 2k ~ 2 \9i e [0, 2tt)1 <j<h—l,8j€ [0, tt/2]}, where we defined as follows: 



COS 9i 

Ik 

e lV k + 1 sin Q i sin 2 CQS 3 



\ 



e Mk sin 9 1 cos 9 2 



sinful sm 02 sin #3 • • • sm%_2 cos 9k-i 
y e i»2fc-2 s i n sin 6*2 sin0 3 • • • sin6> fc _ 2 sm^ fe _! / 



(6) 



The invariant measures v{dff) described above is from [21, p. 31]. 



v{d9) 



r fe-l 



■ Sill 



isin 2/c 5 2 • • • sin0fc_i cos#i cos# 2 • • -cosQk-i d6\ d9 2 ■ ■ ■ d9 2 k-2- (7) 



Lemma 1 If the deviation measure W is characterized as W(p, p) = ho df s (p, p), we can describe 
the maximum possible error of the optical measurement Tl n as: 

V w ^ (n„) = 2(fc - 1) f " + k _ ~ ^ J 2 h(9) cos 2 " +1 9 sin 2 ^ 3 9 d9. 

For a proof, see Appendix [b]. 

Next, we asymptotically calculate the error of the optimal measurements ±I„ in the third order. 



G 



Theorem 4 When the deviation measure W is described as W — d^ , we can asymptotically cal- 
culate the maximum possible error of the optimal measurement as: 

lim = (8) 

Specially in the case 0/7 = 2, we have 

V d2 ^ n \n n )n=^-^ = (k-l)jr(-^y ^k-1 as 00. (9) 



i=0 



When the deviation measure is defined by the square of the Fubini-Study distance, we can asymp- 
totically calculate the maximum possible error of the optimal measurement as: 

P4.'(«)(n„)n S (k - 1) - -k(k - I)- + k(k - l) m ~ 7 i- as 00. (10) 

3 71 45 

TVie error of the sequence {H n }^ =1 of the optimal measurements can be calculated in the sense of 
large deviation as: 

-log (PiH ) {peV(n)\dfs(p,p)>e} 



n 



= logcos^ e+(k- 2)— ^— + (- log(fc - 2)! + 2(fc - 2) log(sine) - 2 log(cose)) - 
n n 

/ k 2 — k — 2 \ 1 

+ h (fe - 2) cot 2 e -r asn-KX), (11) 

\ 2 / 

where Prf^ £? denotes the probability of B with respect to the probability measure tv{M{dui)S) for 
a Borel B cfl, a measurement M € M.(Cl,T~C') and a state S € S(Tl'). 

For a proof, see Appendix The first term of the right hand of (|ll| ) coincide with the logarithm 
of the fidelity. About the fidelity, see Jozsa p2| . 

In this paper, e in equations (|li"| ) is called admissible radius. 

Since 

lim^^ = -l, (12) 
we obtain the following large deviation approximation. 

lim lim J-log(pr^{pe V(H)\d fs (p, p) > e}) = -1. (13) 

e— »0 n— >oo e n V " / 

5 Semiclassical measurement 

In this section, we consider measurements allowed the quantum correlation between finite samples 
only. A measurement M on J{( nm ) is called m-semiclassical if there exists a estimator T on the 
probability space V(H) x • ■ ■ xV(H) whose range is V(H) such that 



M{B)= U m {dp l )®---®n m {dp n )VB(zV{H). (14) 

JT- 1 (B) s v ' 

n 

We compare the error between m-semiclassical measurements and the optimal measurement IT nm 
for nm samples of the unknown state as the equations (p^|) ,(11), (|l3|) . 

In doing this comparison, we bear in mind asymptotic estimation theory in classical statics. In 
classical statics, it is assumed that the sequence of estimators satisfies the consistency. 
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Definition 4 A sequence 

| T (n)|oo =i $ estimators on a probability space Q is called consistent with 
respect to a family {pe\0 £ 0} of probability distributions on if it satisfies the condition 
where every 

Tin) 

is a probability variable on the probability space Q X • • • X Q whose range is 9. 



^ ] {dj{T {n \6) > e} -> as n -> oo,V6> e 0,Ve > 0, (15) 

n ■ (n) 

where dj denotes the geodesic distance defined by the Fisher Information metric and p g denotes 
the probability measure pg x • • • x pe on the probability space O x • • • x O. 



It is well known that the following theorem establishes under the preceding consistency |23|, |2J, £5| . 

Theorem 5 If a sequence {T^}'^L 1 of estimators is a consistent estimator with respect to a family 
{pe\0 S 0} of probability distributions on a probability space f2 which satisfies some regularity, then 
we have the following inequalities. 

lim n J J ■■■ J d 2 j{r^{xi,x 2 ,... ,x n ),e^pf\dx x ,... ,dx n ) > dim© (16) 

n 

lim -log(p^{D(p T(n) \\p e )>e}) > -e (17) 

lim lim ^- log (p^{dj(T( n \ 6) >e}) > -\, (18) 

where D{p\\q) denotes the information divergence of a probability distribution q with respect to 
another probability distribution p defined by: 

D (p\W)-= / (logp(u)-logq(u)))p(u))dw. 
J a 

Under some regularity conditions, the lower bounds of ) can be attained by the maximum 

likelihood estimator (MLE) . A regularity condition for attaining (Y^) is different from a one for ([7 



The lower bound of can be attained by the MLE when the family {pe\0 € 0} is exponential, 
but generally cannot be attained. 

For the comparison, we apply Theorem [|to the family of distributions {trII m (c? / o)p(' m '|p S V(H)} 
given by the measurement H m and the family of states {p^\p G V(7i)}. Let TL m ) be the 
measurement on H'™' defined by the estimator TW and n data given by the measurement 
n m ® • • • ® IT m and the state p( nm \ We consider the sequence of measurements {T^ n m ^\^Li. From 

n 

the symmetry of'P(W) and H m , the information divergence of a probability measure trII m ((iyo)/9^ m ^ 
with respect to another a probability measure tr Tl m ( dp)p^ is determined by the the Fubini-Study 
distance e between p\ and Thus, the divergence can be denoted by D m (e). From Lemma |^, 
the geodesic distance dn m with respect to Fisher information metric in the family of distributions 
{trn m ( dp) p^\p € V(H)} is given by: 

du m = \/2mdf s . 
Since dimViH) = 2(k — 1), we have the following inequalities: 



lim nmT> d f"( nm \T (nm) ) = lim max nm d 2 f Ap,p)tr{T (nm) {dp)p {nm) ) > k - 1 (19) 

n-»oc ' n->oo p£V[H) Jv{U) 

1 , „ „(»"•) , . £> " , Ml 



lim —fog P4 ,{peP(W)|d/.(p^)>e} > ^W) 

n—>oo nm i»>"») m 

lim i Jim * l og p r f"j {p eV(H)\d fs (p,p) > e} > -1. (21) 
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The lower bound of ([l9|) is consistent with the first term of the right hand of ( |Io|) and the lower 
bound of ( pl|) is consistent with the right hand of ( 13 ) . The family of distributions {tr H m ( dp) p^\p € 



V(H)} satisfies a regularity condition for (|16|) by MLE. But, we cannot show that it does a one 
for (jlj). We have the following lemma concerning the comparison of the lower bound — D "^f > of 
(p3) and the first term 2 log cos e of the right hand of dlT 



Lemma 2 We can calculate the divergence D m (e) and the distance djj m as: 

-lAn(e) Sm e I 2 \ 2 

= > — ; >— logll — sin e) = — log cos e as m —* oo (22) 



m 

i—i 



d]j m = \ // 2mdf s . (23) 

Therefore, Dr ^ € ^ is monotone increasing with respect to m. 
For a proof, see Appendix (22) derives that 



0< ' mlOgCOs2 2 e " A " (£) ^0as^0 (24) 
me 2m v ' 

which means that the first term of (|ll]) cannot be attained by a semi-classical measurement. 
However, it is an open problem as to whether the left-hand side of ([l3|) can be asymptotically 
attained by a 1-semiclassical measurement. Concerning MSE, the first term of ( |l0| ) can be asymp- 
totically attained by it i.e. it can be asymptotically attained by measurements without using 
quantum correlations between samples. Thus, in order to attain it asymptotically, it is sufficient 
to physically realize the optimal measurement IT on a single sample. Indeed, IT is a random 
measurement as follows. To denote IT as a random measurement, we will define the simple mea- 
surement E g (g S SU(fc)) whose measurable space V(H). For an element g e SU(fc), the vectors 
4>i{g), ■ ■ • ,4>k{g) in Ti. are defined as: 



(<Mff) • • • <MsO) =9- 



The measurement E g is defined as: 



E g {\U9)){H9)\) = \Hg)){H9)Y 

Therefore, the optimal measurement IT for a single sample can be described as the following 
random measurement: 



IT = / E g p{dg), (25) 

JSV(k) 

where p, is the invariant measure on SU(fc) with /i(SU(fc)) = 1. Therefore, in order to realize the 
optimal measurement IT, it is sufficient to realize the simple measurement E g for any g € SU(/c). 



6 Conclusion 

We have compared two cases. One regards the system consisting of enough samples as the single 
system, the other regards it as separate systems. Under this comparison, the MSEs of both cases 
asymptotically coincide in the first order with respect to the Fubini-Study distance (see (|l0|) and 
(|l9|)). However we leave the question of whether they asymptotically coincide in the second order 
with respect to the Fubini-Study distance to a future study. On the other hand, in view of the 
evaluation of large deviation, if the allowable radius is finite, neither coincide (see ( [ll] ) and (p0|)). 
However, in the case of the allowable radius goes to infinitesimal, it is an open problem as to 
whether both coincide (see ( |l3| ) and (|2~l|)). 

These results depend on the effect of a pure state. Therefore, it is an open question as to 
whether the MSEs of both cases asymptotically coincide in the first order in another family. In 
the case of large deviation, the same question is also open in the limit where the radius of allowing 
error goes to infinitesimal. 
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Appendices 



A Proof of Theorem 3 

In this Appendix, assume that p = |0(O))(0(O)|. Because Tli n ^ is irreducible with respect to the 
action of SU(fc), 



= (e^° 



(n) 



Oi € C,gi e SU(fc) 



(26) 



We assume that W{p,p) = h(trpp). As h is monotone decreasing, there exists a measure h! on 
[0, 1] such that h(x) = h'([x, 1]). 

The function hp on [0, 1] and the deviation measure Wp are defined as follows: 



hp(x) 



1 for x < (3 
for x > (3 

Wp(p,p) := hp(trpp). 
From Lemma [|, for any measurement II we have 

V^ n \U)= [ vJ^ {n) (U)h'(dp). 

J [0,1] 

From (|26|), it is sufficient to show the following for {(f>i\ C TL in the case of W = Wp. 



tr Wp(p) 













From Lemma |] it is sufficient for ( |27j ) to prove the following: 



E 



E^V^"' Id-W^Uo)^ 



> (m in) \wp(p)\m in) ) ■ (Y,^ 



O) 



Id-Wp(p) 



E 



Remark that | (0(0)10(0)) | 2 = cos 2 0i. From Lemma ||, we get 



E 



Wp(p) 



E 



/,(») 



Id-W^p) 



E^ 



fc'- (fc-1)! p 



^(fe-i) 

fc ; ■ (fc-i)! r 

TT^- 1 ) Jo 



/i (0i ) cos 0i sin 2 *" 3 0i d0i 
/i (0i ) cos 0i sin 2fe - 3 0id0i 



0(O) (n) W-^(p)U(0) w ) = C/ cos 2n+i 0i sin 2fe - a 0i d0! 



(«)\ _ 



0(O) (n) Id-W^p) 0(O) (n) 



C / cos 2n+1 0isin 2fc - 3 0id0i, 



(27) 



(28) 
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where 



,6k-i) 
C 




h,... ,e k _ 1 )\{d6 2 --- de k _ x ) 
Y,(<f>i\<f>(0)) n (<t>{0)\<t>j) n de k --- d9 2k . 



1,3 



fc-1 



fc'-(fc-l)! 



2lX /.27T r% 




X(d9 2 d9 3 --- dS k -i)d0 k --- d9 2k - 



X(d9 2 ,... , d6 k -i) := sm 2k 5 6 2 ■ ■ ■ sin6» fc -i cos6> 2 ■ ■ • cos6> fe _i d6 2 ■ ■ ■ dO k -i. 



Therefore, it is sufficient for the equation (Eq) to show that for tt/2 > 9\ > 9[ > 



/i(0i) sm^ 6»i cos Zn+l 9[ sm^- A 9[ > f 1 (9[)sm 2k - 3 9[ cos 2n+1 V x sm~ r/, . 
It suffices to verify that for 9i G [0, f ], 2 < i < k - 1, tt/2 > 6>i > 9[ > 

h{9 u 9 2 ,... ,d„_i) > f2(9[,9 2 ,... ,9 n - x ) 



2k-3 , 



COS 2 " 6»i 



COS 



2n fl/ 



Thus, it is sufficient to prove that the following is monotone decreasing about 9\ for any 9 2 , 



1 



COS 2 ™ 9\ J 



J2(&\m) n (m\h) n d9 k ■ ■ ■ d9 2k _ 2 



k-l 



Letting 



"Ptrh 2 



we get 



cos n 6»i 



k-l 



-(pi + e^ 6 "-^-^ tan 9 X sin6> 2 • • • sin^ cos9 j+1 (, 



i=2 



+e i(9 2k -2-^ z ) t an6»isi 



smt/ 2 ■ ■ ■ smt>k-i(t 



■>k-X- 



(29) 



Letting x :— tan6*i, Lemma || induce that ( |29| ) is monotone decreasing about 9\. The proof is 
complete. 



Lemma 3 // the deviation measure W{p,p) = h'([tr pp, 1]), i/ien 

y[o,i] 



(30) 



ii 



Proof For the probability measure 7r on V(Tl), we have 



W(p,p)ir(dp) 



V(H) 



h{tr Pp)tt( dp) 



V{H) 



V{H) J [0,1] 



hp(ti pp)ti(d/3)n{dp) 



[0,1] ^JV(H) 



hp (tr pp)Tr{dp))h'(d/3) 



1 (/ Wf,(p,pMdpj)h'{d0). 

[0,1] K JV(H) ' 



Substituting Tr(dp) for tr(II( dp)p^), then we obtain (|3(]|). 



□ 



Lemma 4 Le< 7i 6e any finite dimensional Hilbert space. For any elements tf>, if) € "H and any 
selfadjoint operator A on Ji, the following are equivalent. 



Lemma 5 we /iawe 



E 



1^3 (j>) 



E 



E^ n) 
* / 

E^ 



*' • (Ar-l)! 
fc' ■ (fc-1) ! 



/i (^i ) cos 0i sin 2fc ~ 3 0id0i 
A (0i ) cos 0i sin 2A: - 3 0i d0i 



(<£(0) (n) Wis(p) U(0) (n) ) = C/ cos 2 ™ +1 0i sin 2fe - 3 0i d0i 
^(O) (n) |ld-W>(p)|0(O) ( "^ = C cos 2n+1 0isin 2fc - 3 0id0i. 
Proof W^(p) is denoted as follows: 



JV{H) 

= k' [ P (n) v(dp). 

J {p£V(H)\tr pp<(3} 



We obtain 



W (p) 



E 



fc' / P {n) v{dp) 

{p£V(H)\trpp<l3} 



{pev(H)\u pp<0} 



(n) 



E 



(n) 



{pe-PCH)|trpp</3} 



fc'-(fc-l)! 
Tr(fe-i) 



/i (0i ) cos 0i sin 2 ' 1 " 3 0i dBi 
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Similarly, 



E 



Id -Wp(p) 



E 



(f){0) in) ld-W fj (p) ^(0) ( 



1,3 



~~" <{p£V(H)\trpp>f3} 



k' ■ (k- 1) 



7T 



(fc-1) 



: / /i(Oi) 
Jo 



cos6»i sin 2fc - 3 0i d6»i 



(0(O)W W (p) m {n) ) = k' f 
x 1 ho 



k' i (4>(o)\p\m) n Hdp) 

{p€V(H)\tvpp</3} 



= C / cos 2n+1 6>i sin 2 *" 3 1 d6x 



k' / (<f>(0)\p\m) n v(dp) 

J{peV(H)\tvpp>f3} 

C [ cos 2 ™ +1 6» 1 sin 2fe - 3 l( i6' 1 . 
Jo 



Lemma 6 The following function f(x) is monotone decreasing on [0, oo): 



: =EE r ■■■ r ( c ^ d ° a + x E n Ue< + * E e-^ci) 1 

a=16=l^ ° ■ 3 = 1 3 = 1 



fix) 



where c 3 n ,d J n are any real numbers. 
Proof The set K™ is defined as follows: 



1 1 =(/!,■■■ ,/ m )e(N+'°r 

The number C(J) is defined for I e X™ as sufficing the following condition: 



Therefore, 



cg + ^c^^l = J] C(/) e id "( C 0) /o e l/l ( ei+d ^( C i) /l ... e ^^+' i ")(^) / ^ n - / ° 
i= Y I ieK* 



r k+i 



Thus, 



= EE^E^) 6 ^"^^ 

a=-l b--l I 

m m 

= (2n) k J2 C(I) E E eK12U ^4) (c O)/o . . . {c ky k {c oy„ {c ky kx 2n-2i 

I a=l 6=1 

= (2 7 r) fe ^C(/)^(/)x 2 



a=l 6=1 

r 2n-2/ 
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where 

m m 

m ■= E E ^ u h ^ u Udl) ( c ') /0 • ■ ■ i<) Ik ■ ■ ■ (c k b ) Ik ■ 
0=1 6=1 

It is sufficient to show D(I) > 0. Letting 

Va : = (c o a y°...(c k a y* 

k 

y a ■■= E 7 ^ 

w a ,b ■■= cos(y Q - y b ), 

we have 

rn m 

D(I) = E E V a w a,bVb- 
a=l b=l 

Then 

w a<b = cos(y Q - yb) = cosy a cosy b + sin y a sin y b . 

As {cosy a cosj/f,} and {sin?/ a sinj/b} are nonnegative, {w a ,b} is nonnegative matrix. Therefore, we 
obtain D{I) > 0. □ 



B Proof of Lemma [j] 



V w >W(IL n ) = / &(d /a (p,p)tr(II„(dp)pW) 



■P(H) 



( U t k , 1 I cos 2 " +1 0! sin 2fe - 3 0! d& 

1 fc — 1 / 7r K-i 




sin 2fe 5 6*2 • • • sin6» fc _i cos6> 2 • ■ ■ cos6» fe _i d6> 2 • • • dd 2 k~2 



11 



fe-l fe-2 



\ fc — 1 / 7T K_1 

x / x 2k ~ 5 dx ■ ■ ■ [ xdx-^nf- 1 
Jo Jo 

= 2(k-l)( n+ k k _~ 1 ^j Jj h(8) cos 2 ™ +1 sin 2fc - 3 d0. 



The proof is complete. 
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C Proof of Thereon § 

Definition [l] and Lemma [l] means that: 

Ztf.M (n „) = 2{k - 1) ^ + k _ ~ J 7 cos 2n+1 sin 2fe - 3+ ^ d0. 



Since 



we have 







cos x 6> sin y 6» d<9 = v 2 ^ v 2 ; Vs. y € R. 



n + fc - 1\ T(n + l)T(k - 1 + 7/2) 



^*>(n„) - ^"'l. , r(n + t + 7 /2 l 

r(n + l)r(fc - 1 + 7/2)r(n + fc) 
r(n + fc + 7/2)r(n + l)F(fc - 1) 

r(n + fc) r(fc- 1 + 7/2) 
r(n + fc + 7 /2) r(fc-i) ■ 

Therefore we obtain (|J) from the following formula of T function: 

lim F( " + X) = 1. 

n— >oc T(n)n x 

Letting 7 :— 2, we obtain 

r&Mm 1 = T(n + fc) F(fc-1 + 1) = k-l_ 
K n> T(n + fc + l) r(jfe-l) n + fc' 

Thus, we get (0). 



(2i - 2)!! sin 2i 6> 



Next, we will prove (|T^). 9 2 can be expanded as 

00 

* 2 = E 



-(=0 



(2*-l)!! i 



where we put (2n)H = 2ra(2ra - 2) • • • 4 • 2, (2n - 1)!! = (2n - l)(2n - 3) • • • 3 • 1, 0!! 
From (|32l ) we have 



P4 s .(") (n „) = g^^|^r.(«) (n „) 



i=0 

i-1 



^ (2z - 2)!! tt * - 1 



(2z - l)!!i 11 n + k + j 

i=o y ' j=o J 

- (k- i)l -h(k- i)\ + k {k-i)^=l±. 

n 6 n z 45 n a 

Thus we obtain ([n|) . Lemma [l] derives that 



l°gPr n „ {p&V(H)\d fs (p,p)>e} 

k - l x '" cos2 e 
fc- 1 



log I (fc- 1)( "7" M / x n (l-x) k ~ 2 dx 
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Therefore, it is sufficient for ( |l3| ) to show that 
log 



n + k — 1 
k - 1 



(k - 1) logn - logffc - 1)! + — ^ 

n 2 



log 



(/>cos 2 e \ 
J x n (l - x) k ' 2 dx J = 2nlogcose - logn + 2(fc - 2) logsine - 21ogcose 



--(l + (fc-2)cot 2 e) 
n 



The left hand of (fM) is calculated as: 



log 



k-l 



n + k - 1\ 

7 i=0 



n + i 



k-l 



= (k-l) logn - log(fc - 1)! + V log(l + -) 

* » T? 



fc-1 . 

2 



(fc-l)logn-log(fc-l)! + V- 



(fc-l)logn-log(fe-l)! + 



1 (k-l)k 
n 2 ' 



Therefore, we have (34). The left hand of (35) is calculated as: 



/ /"COS 2 € \ 

gij x n (l - x) k ~ 2 dx \ - 2nlogcose 
Vcos^e/ / 



lo; 
log 

log ( / x n (1 



2 \fe- 2 1 j 
cos a; I — dx 



'k-2 



= —2 log cos e + log ( | . ^ ) ( — cos 2 e^ 



n + i + 1 



2 log cos e- log n + log f ^ ^ . ^ (-cos 2 e)* j^j+i 



\i=0 
'k-2 



-2 log cose — logn + log I ^ 



k- 2 



2 \* 



cos e 



) I" 



i + 1 



•2Iu,!i<-osi 1< >:;// ! 1<>:; j (l - cos 2 e)* 2 - - (l - cos 2 e) k 3 (l - (k - 1) cos 2 e] 

2 \ k ~ 2 i A 1 1 - (A: - 1) cos 2 e 

= — 21ogcose — logn + log (1 — cos e) + log 1 

x ' \ n 1 — cos z e 

i 9 \ 11 — (k — 1) cos 2 e 

= -2 log cose - logn + (k - 2) log (1 - cos z e) . 

v ' n 1 — cos^ e 



We obtain <M). 
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D Proof of Lemma 

From the symmetry of Vi^H) and n m , we may assume that p\ = |</>o){</>o|i P2 = \4>e){4>t\- First, we 
consider the case of fc = 2. For the following calculation, we prepare the following equations: 

m\m)\ 2n 

= (cos 2 e cos 2 6\ + sin 2 9% sin 2 e + 2 cos e sin e cos 0i sin 6% cos 2 ) n (36) 

/ log(l + 2acos6» + a 2 )d6» = 47rVj(|a|)log|a|, (37) 
Jo 



where the function ip is defined as: 

ip(x) = 

Paying attention to (g) and (]?]), we have 

— -Dn m (px ^11/4"^) — to log cos 2 e 



1 z > 1 

a; < 



1 , ,,J , l(0ol^W)l 2m 
-l (m + 1) i(«) l0g ra(^ 



/ log V I m 2m K^ok(g))l Kdg) +mlogcos 



2(m + l) 



7T 

log 



2 ^ r * ' cos 2 0i cos 2 e 



Jo 



(cos 6*i cos e + sin 9\ cos 02 sin e) 2 + (sin 0i sin 62 sin 



■cos 2m+1 0! sin0! d9 1 d9 2 
(m + 1) 



7T 

p27T 



/ 2 / log (1 + 2 tan 0x tan e cos 2 + (tan 0! tan e) 2 ) d0 2 cos 2m+1 0! sin0! 
Jo Jo 

/ 47T log(tan 0i tan e)-0(tan 0i tan e) cos 2m+1 0i sin 0i d0i 
Jo 

J 4tt log(tan 0i tan e) cos 2m+1 0i sin 0i d0i 



/o Jo 

(m + i) n 

* Jo 

(m + i) n 



= 2{m+\) I log(tan 2 0itan 2 e)cos 2m+1 0isin0id0i 



= (m + l)J log -^tan 2 ej x" 1 ^. (38) 

Substitute a = tan 2 e in (^l|) of Lemma || , then 

/■sin 2 e / 1 - x \ m sin 21 e 

(to + 1) / x m log tan 2 e dx = - log cos 2 e - V — : — . (39) 

Jo V x J i=l 1 

From (|38| ) and ([39]), we have 

£>n m U m) ||^ m) ) _^sin 2 * e 



E— • (4°) 



Therefore, we can prove (E2 
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Next, we consider the case of k > 3. In this case, we have: 

m\m)\ 2n 

= (cos 2 e cos 2 9\ + sin 2 0\ cos 2 9 2 sin 2 e + 2 cos e sin e cos 9\ sin0i cos 02 cos0fc)™. (41) 
Paying attention to (|J),(]37|) and Lemma 0, we can calculate as: 

-Dn m (p[ m) \\p { 2 m) ) -™ log cos 2 e 
m 

- 4 (C"I- r i<*i*»r**> + *.*~«) 

2(fc-l)(fc-2) (m + k-l 



TV 



k - 1 



27r rZ / cos 2 ^ 1 cog 2 £ 

log 



JO Jo 



cos 2 e cos 2 0i + sin 0i cos 2 6*2 sin" e + 2 cos e sin e cos 9\ sin 0i cos 6*2 cos 9k ) 

k 



• cos 2m+1 6>i sin 2fe - 3 6»i cos 2 sin 2fc " 5 9 2 d6 x d9 2 d9 



2(fc- l)(Jfe-2) (m + k-l 



7T 

f-27T 



(m + k — X\ 



o Jo 



■ / log (l + 2 tan 0i cos 2 tan e cos 6*^ + (tan0i cos 02 tan e) 2 ) d9k 
Jo 

• cos 2m+1 0i sin 2fc_3 0i cos 2 sin 2fe_s 2 d9 x d9 2 
2(k- l)(fc-2) (m + k- x\ n n A ,, 

= -i ^ M k 1 ) J J 4?n/>(tan X cos 2 tan e) 

• log (tan 0i cos 2 tan e) cos 2m+1 9 1 sin 2k ~ 3 6i cos 2 sin 2 ' 1 " 5 2 d6 ± d9 2 . (42) 

Substitute that a := tan 2 e, s := sin 2 2 , y := cos 2 0i, then the condition 
tan 0i cos 02 tane > 1 turns into the following conditions: 

ii a 

1 - " s > x > 0, — — > y > 
a(l — y) 1 + a 

Using (fi"2| ) and ((5^), we have 
— L)u m [pi \\p 2 ) — to log cos e 

TO 

= (*-!)(*- 2)("' + * " ') (f^ x'"> log ((1 - *)^^) <k) iTd - dy 

m + k-l\ f~ , ( y \ m/H Nfe _ 2 , ^ (m + k - 1 



' log — ^h/ m (l-2/) fc - 2 dy-(fc-l) ' 1 / 



fe — 1 / 7 ya(l — y)/ \ fc— 1 / \l + a 

(43) 



where /(#) is defined as: 

fe-2 
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From Lemma g, the derivative of f(x) can be calculated as: 



'k — 2 / y \ i — 1 ^ 



\ 1-2/ 



/'(*) = y n ^ (E(i-f ) J» m (i-y)*- 3 * 



\ 



1 - 



1-g 



1 


ar(l - 


X) 


1 




x(l - 


x) 


1 




x(l - 


x) 






x m 




x{l- 


X) 



y m (l-y) k - 2 dy 



y m (l-y) k - 2 dy 



x(l — x) J Q \ X 



y 



fc-2 



i-- y md y 



fc-2 

E 

i=0 
fc-2 

E 

i=0 



fc-2 



-l rl 



(-W+'dy 



(1 - t) fc " 2 * m 



fc - 2\ (-a;) 4 x 



m+l 



i J m + i + 1 x(l — x) \ fc— 2 / m+l 



fc-lV 1 1 



By ( |5l|) and Lemma [| the first term of (43) is calculated as: 



-(fc-1) 
= ~{k-l) 

= "(fc-1) 



m + fc — 1\ /" 1+ 
fc-1 



y m {1 _ y) k-2 log 



y 



a(l-y) 



dy 



fc — 2 

m + fc — 1 \ v— v /fc — 2 



i=0 
fc-2 



2/ 



m + k- 1\ ++ (k-2\ (-1)' 
fc-1 / frf 



i=0 
fc-2 



0(1 "») 

m+i 



i J m + % 



T l-log(l + a ) + £i( T a 



m + fc- 1\ ^ ^fc-2^ (-I) 1 
fc-1 7 ^ 



i=0 

1 / 'I 

J 



i I m + i + 1 



j m-\-i 



1 / a 



j=l J 7 j=m+l 

fc-2 



; )(V 



fc- 



fc- 1 



fc 



= "(fc"!) 

-(fc-1) 

= "(fc"!) 



si=0 
fc-2 ', 



fc-2\ (-1)' \ / , n . ^ 1 / a 



i EE 



fc-2 



(-1)* 



j'=i 



* / (m + i + l)(j' + m) \1 + ' 
m + fc — 1\ f /m + fc — 1\ 1 1 



fc - 1 



m + k — 1 



fc-2 



1 / a 



T -log(l + a) + £- I 



fc-1 7 V 1 + a 



1 / a 



3 = 1 



(fc-1) 



m + fc — 1 



fc-1 /Vl+a/ 
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: =EE 



k - 2\ 



where g(x) is defined as: 

fc-2 i 

— ' ^— ' V i / (m + £ + llfi + m) 

2 — J — 1 V ' 

By Lemma ^, the derivative of g{x) is calculated as: 

fc-2 



fc 2\ ( jT \ " r ; /.« J 

'' 2 'k-2\ (-1)* _ x . m l-^ 



i / (m + £ + 1) f — .t 



E 

i=0 

£ /(m + i+1) l-x2-*>\ i J (m + i + 1) 

x m /to + fc - l\ _1 f a;" 1 y?/fc-2\ (-a?)* 

1 - x \ fc-2 / to+1~1-x^V i / (m + i + 1)' ' ' 



Prom (Q) and (ff6|), we have /'(a;) = —g'(x). The definitions of f(x) and 5(2;) means that /(0) 
g(0) = 0. Then we obtain /(a;) = -£f(x). By @ and (§|), we have 

-^n m (Pi m) ||p 2 m) )-2TOlogcose 



/, ^fm + k — 1\ / ( a 
-(fc- 1) 



- 1 A U+fl/ V 1 + a 



1 / a 



= 1 °g( 1+a )-E-( v T +a 

= — log cos 2 e — - 
Then, we obtain 



sin 2 - 7 e. 



1=1 J 



n / (m)n <m)\ m . 2i 
£> n,„(pj lPg ) = sm e 

to i 
3=1 



We proved (|22|). 

Next we will prove (23). We consider the tangent space TpV(H) at p := \cf>(0))(<f)(0)\. If c(t) 
is a curve on V(H) such that c(0) — p, c denotes the element of T P V(H) defined by c(t). the 
Fubini-Study metric gf s is defined as: 

,. , /,. d /s (c(0),c(i)) N 2 
Therefore, it is sufficient to show that 

J n n = 2n 9fs- 
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Let c(t) :— \(p t )((f>t\, <t>t '■= <f>(t, 0, . . . ,0). (See the equation (||).) Because g/ s (c, c) = 1, it is sufficient 
to prove that 

(c, c) = 2n. 

We assume that fc > 3. From (^d]), we have 



■^log(|(0 t |^))| 2 ")' 



t=0 



K^ol^^))! 2 " = 4n 2 cos 2n - 2 0isin 2 0icos 2 02cos 2 fe . (47) 



By @ and @, we have: 



m + fc — l\ /" / d 



log(|(^|^))| 2m )L =0 K0o|0W)| 2m i/(df) 



2(ft-l)(fc-2) (m + k-1 



/m + fc - 1\ 9 /" 



4m 2 / cos 2 " -1 6i sin 2 *- 1 X d0 x 



2tt 

cos 3 02 sin 2fe ~ 5 2 d0 2 / cos 2 9 k d0 k 
o 



2(fc- l)(fc- 2) fm + fc- l\ 2 (m - !)!(&- 1)! l!(fc-3)! 



7T V / 2(m + fc-l)! 2(fc-l)! 

= 2m. (48) 

We get ([23"1). In the case of fc = 2, similarly we can prove (p3j). 

Lemma 7 If k > 3, then we have 

f(d u e 2 ,e k )u(de) 

V(H) 

= 2 ( fc - 1 )( fc ~ 2 ) f f 2 f 2 f(0 1: e 27 9 k )cose 1 sm 2k ~ 3 e 1 de 1 cos9 2 sm 2k - 5 e 2 d9 2 de k {49) 
k Jo Jo Jo 

Proof From (Q) the left hand of ( fl9| ) is calculated as: 

f(ei,e 2 ,e k ) v (de) 

V(H) 

^~^ ! / / 2 I* f(9 1 ,9 2 ,9 k )cos9 1 sin 2k - 3 9 1 d9 1 cos9 2 sm 2k - 5 9 2 d9 2 d9 k 

K k Jo Jo Jo 

f rf / ,27r /" 27r 

• • • / sin 2fe " 7 3 • • • sin fe _i cos 3 • ■ ■ cos fc _i d0 2 • • • d0 fe _i / • • • / d0/c+i • ■ ■ d9 2k _ 2 
Jo Jo Jo 



k-3 k-1 

2tt 



2 / 2 



/(0i, 02, Ok) cos 0i sin^ 3 0i d0i cos 2 sin^- & 2 d9 2 d9 k 
it Jo Jo 

(fc-1)! f ju_ 7 / , , _\fc-2 



ji- / x^-'das--- / jcda;-(27r)* 



fc-3 

27r rf rf 

/ / f{9 1 ,9 2 ,9 k )cos9 1 sm 2k - 3 9 1 d9 1 cos9 2 sm 2k - 5 9 2 d9 2 d9 k 
Jo Jo 

(fc-1)! (27r) fc - 2 

x n k-i 2 fc -3(fc-3)! 

2(fc-l)(fc-2) ^ 



7T 



/ / 2 / 2 /(^I^2,0fc)cos0isin 2fc - 3 0id0icos02sin 2fe - 5 02 d0 2 d0fc. 
Jo Jo Jo 



Then we obtain (Ea). □ 
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Lemma 8 The following integral can be calculated as: 



pl-a / 1 _ \ 1 / m+1 

i^^fe))^ ' STl(- l0 « (1 + « , + §KiTi)')- (51) 



Proof The equation (|50| ) is derived by the following: 



a / m+i , n 

/ x m log(l — — (a m+1 - 1) log(l - a) - V - 



(52) 



Also, the equation ( |5l| ) is derived by ( |52| ) and the following: 



x m log a; rfx = — !— ( a m+1 ( log a ■ ) ) . ! .18 ) 

m + 1 \ V to + 1 



□ 



Lemma 9 VFe /lave i/ie following equations: 



±( n Y-^= [\™-\i- x rd X = ( m+n Y± (54) 

^ \ i ) m + i J v ; \ n J m 



i=0 

It is easily derived. 
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